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Abstract 

The structure of ^Li, ^Li and ^^Li nuclei is investigated in a model space which 
includes all configurations with oscillator energy up to 3huj above the ground 
state configuration. The energy spectra and electromagnetic properties of the 
low-lying states are determined with various two-body interactions, which are 
derived from the Bonn potential. The results of these shell-model calculations 
depend on the strength of the tensor component contained in the NN interaction 
and also on the treatment of the Dirac spinors for the nucleons in the nuclear 
medium. In addition the calculation determines the dipole polarizability of ^Li 
and ^^Li caused from virtual El excitations to the positive parity states of these 
nuclei. It is demonstrated that the BAGEL scheme provides a very powerful 
tool to consider contributions from virtual excitations up to large energies. 
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1 Introduction 



The static and dynamic electric moments of nuclei are traditionally measured by 
Coulomb excitation^). The quadrupole moments of the excited nuclear states can 
be measured by the reorientation method^). It has long been established that the 
quadrupole moments as well as the B{E2) values can be affected from virtual exci- 
tations to the giant dipole state^). In particular, if the dipole state occurs at low- 
excitation energy, it can cause a measurable renormalization of the E2 operator. 

The most careful measurements of the correction due to the dynamic dipole po- 
larizability on the B{E2) and the quadrupole moments has been in the case of the 
nucleus ^Li ^~^). This nucleus has a ground state with J'^ — ?,/2~ and a first excited 
state with J'^ — l/2~. The excited state has no spectroscopic quadrupole moment 
and hence a precise determination of the B{E2) f value for the excitation of the l/2~ 
state over a range of energies will enable one to measure the inelastic dipole polar- 
izability. Similarly, a careful measurement of the ground state quadrupole moment 
will enable one to measure the elastic dipole polarizability. The results of various 
measurements are reviewed by Barker et al ^) and Voelk and Fick^). 

There exist several theoretical investigations on the dipole polarization effects 
of ^Li. Several of these calculations adopt the cluster-model approach^"^^). On 
the other hand the shell-model approach has been followed by Gomez-Comacho and 
Nagarajan^^) and Barker and Woods^^). In both these calculations the natural parity 
states are determined considering only Op shell configurations while the unnatural 
parity states are studied in the space of the Ihcu excitations. 

Recent experiments^^) suggest that a component of the dipole resonance, termed 
the soft dipole mode'-^) or the pygmy resonance'-^), occurs at a very low excitation in 
^^Li. This feature combined with the large r^s radius of 3.16±0.11/m ^^) have made 
this nucleus the subject of several theoretical investigations^^"^^). In these studies 
various approaches have been adopted, hke the cluster- model^^"^^), the random phase 
approximation^^), the sum-rule approach^^) and the mean-field description^^"^^). In 
addition the structure of ^^Li has been investigated in two recent large scale shell- 
model calculations by Hoshino et al ^^) and Hayes and Strottman^'^) . In both these 
calculations, the properties of the natural parity states are described in the space of 
(0 -|- 2)hu! excitations. For the unnatural parity states these authors employ a model 
space, which contains the Ihcu configurations and some selected 3hu! excitations. A 
common feature in these two shell-model calculations is the manner in which the 
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effective hamiltonian is chosen. This consists of a mixture of hamiltonians each of 
which was originally constructed for calculations in smaller model spaces. Thus Hayes 
and Strottman^'^) use the Cohen-Kurath^^) p-shell matrix elements, the Millener- 
Kurath^^) particle-hole interaction and single particle energies and the Chung^°) sd- 
shell matrix elements. A similar approach is also adopted in the calculation of Hoshino 



In this paper we present a shell-model study of the natural parity states of Li, 
^Li and ^^Li. The properties of these states are investigated in the space of (0 + 2)huj 
excitations. For '^Li and ^^Li we also examine the unnatural parity states for which 
we consider the complete space of (1 + 3)huj excitations. The most detailed data exist 
for ^Li and wc use them to draw conclusions about the reliability of our model. On 
the other hand, there exist very little data for ^^Li and thus our results should be 
taken as theoretical predictions. 

The main purpose of our investigation has been to determine the El polarization 
effects on the quadrupole moment of the 3/2^ ground state of both ^Li and ^^Li as 
well as on the B{E2) | value corresponding to the excitation of the 1/2^ state of ''Li. 
Several sets of results have been determined corresponding to different hamiltonians. 
To test the reliability of these calculations we have also determined the properties of 
the low-lying states of ^Li, ^Li and ^^Li and compare the theoretical predictions with 
the existing experimental data. 



of a realistic hamiltonian which contains no adjustable parameters. Thus in our 
calculation we consider harmonic oscillator energies and wavefunctions for the single- 
particle states, while for the two-body part of the hamiltonian we employ G-matrices 
that have been derived from the Bonn potential"^^). Since our model space, although 
large, is restricted, the final results depend on the choice of the harmonic oscillator 
size parameter. We determine an optimal oscillator parameter by requiring that the 
calculated energy for the ground state for the nucleus is a minimum with respect to 
a change of the oscillator basis. Another new feature of our calculation is that we 
take into account a possible change of the Dirac spinor for the nucleons in the nuclear 
medium as compared to the Dirac spinors in the vacuum. Recent investigations on 
open shell nuclei^^'^^) have demonstrated that effects of a realistic spin-orbit splitting 
can only be obtained if this relativistic feature of nuclear many-body theory is taken 



et al 26). 
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into account. Finally, we consider bare electromagnetic operators. We expect the 
renormalization of the operators to arise naturally from the use of the large model 
space. By considering different model spaces and their effect on the predictions of the 
calculation we expect to obtain a better understanding of the relationship between 
our approach and the cluster model^~^^'^^~^^). 

In sect. 2 we discuss the details of our calculation, while in sect. 3 we present our 
results. Finally, sect. 4 contains the conclusions of this work. 

2 Details of the calculation 

2.1 Definitions 

The effective potential which is used to describe E2 excitation of the projectile in 
Coulomb scattering experiments consists of two parts^) : 

= Vcoup + Vpol , (1) 



where 



with 



^Mfki^) = J^^JJ^^2T.y2,{r){J^Ma^^\JfMf) (4) 



In the above eqs (|^)-(^), Z is the charge number of the target, r is the distance 
between the two nuclei, while |/) and \i) denote the final and initial states of the 
projectile, respectively. 

The polarization potential (H) consists of a monopole and a quadrupole term. The 
monopole term is defined^) as: 



9 ^ 
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where 

and \n) are the states which are connected by the El operator with the state 

The quantity Tj/ in the quadrupole part of Vpoi is the tensor moment of the electric 
polarizabihty and is defined'^'^^) as: 



Ti 



The relation between Tj/ and the quantity S{E1) defined by Hausser et al is 



r., = 'iJ^Sm (8) 



Prom the analyses of the ^Li data^'^) values for four quantities have been extracted. 
These are B{E2; 1 —>■ 2), tu, tu and the static quadrupole moment Qs- We use here 
the notation usually adopted^'^) where 1 refers to the ground state and 2 to the first 
excited state. Qg is defined as: 



IGtt 



■^^^^■^^ -{i\\Em ■ (9) 



5 \ (J, + l)(2J, + l)(2J, + 3) 



The experimental values for these four quantities are compared with the predictions 
of the calculation in sect. 3. 

2.2 The shell-model calculation 

In this section we discuss the details of the shell-model calculation on the Li isotopes. 
We consider a hamiltonian of the form 



H = J2U + T.Gij + lAmu;^R''-^huj , (10) 

i i<j 
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where t denotes kinetic energy, G is the two-body interaction while A and R stand 
for the mass number and center-of-mass coordinate of the nucleus, respectively. The 
term ^Amuj'^R'^ has been included in the hamiltonian H in order to remove spurious 
state effects'^^) while the term — f^cu removes the center-of-mass contribution from a 
non-spurious state. The hamiltonian (0) can be written as: 

H = Ho + V , (11) 

where 



^13 ^ mcu^ 

Ho = Y.iU + ^rnuj'rf)-^hu, V = J:iG^, - ^{n - fjf) • (12) 

i i<j 



The basis of our calculation consists of all eigenvectors of Ho which have un- 
perturbed energy not exceeding by 3hu! the energy of the ground state configuration 
(os)^(oj9)^~^. Thus the model space contains 15 single particle orbitals from the Os up 
to the sdg shell of the harmonic oscillator potential. The basis vectors are constructed 
assuming no core state and are represented as 

|$) = |C^;jr/i) , (13) 

where denotes a distribution of the A particles among the single particle-orbitals, 
while the index /i distinguishes the vectors which correspond to the same set of 
{C"^, J,T} quantum numbers. Matrix elements of operators between the states (O) 



have been determined using the generalized fractional parentage formalism of Skouras 
and Kosionides'^''). 

Since the vectors ( |T3|) have been constructed using isospin formalism it is conve- 
nient to consider in the same formalism the operators EX of sect. 2.1. Thus 



where 



A A 



ES = ^T.^-Y,,{h), Ell = -Y^ro{i)rlY,,{h) . (15) 
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Thus the matrix elements of the operator EX can be expressed as: 



{JfTfMrWE^WJiTiMT) 



6TfTAJfTf\\\E^''\\\J^T^) + ( J/Tj 1 1 1 1 1 1 J,T,) (T^MylO |T^Mt) 

{2Tf + 1)1/2 



(16) 



where 

Z - N 

Mt = ^- , (17) 

Z and being the proton and neutron numbers of the nucleus under consideration. 
In (|l^) and elsewhere a triple bar denotes a matrix element which is reduced in both 
spin and isospin spaces. 



As is evident from ([T5|) , the isoscalar part of the El operator is proportional to 
the center-of-mass vector R. Thus the operator connects a non-spurious natural 
parity state to a spurious one, which belongs to the space of the unnatural parity 
states. As a consequence only the isovector part of the El operator is considered 



in determining El matrix elements (16). However, as discussed below, the isoscalar 
part of the El operator is used in the elimination of the spurious unnatural parity 
states. 

As described in sect. 2.1, to determine the matrix elements of Vcoup and Vpoi one 
needs to determine the wave functions |/) and \i) of the final and initial states, re- 
spectively. In addition, as eqs (|) and (^ imply, one needs to determine all states 
|n) which are connected to both \i) and |/) by the El operator. The determination 
of the initial and final states presents no problem. Their wave functions are obtained 
by straightforward diagonalization of the hamiltonian in the space of (0 + 2)huj exci- 
tations. Moreover, since the initial and final states both belong to the low-lying part 
of the spectrum, one needs to determine only few of the eigenvalues, the lowest in 
energy, of the hamiltonian matrix. 

On the other hand there are problems with the determination of the states \n), 
discussed above. The space of their calculation consists of the (1 + 3)huj excitations 
and, consequently, it has a large dimension (over 3000 for some states of ^^Li). In 
addition the relatively few of the states In) which are strongly connected by the El 
operator to the \i) and |/) states need not necessarily be among the lowest in energy. 
Thus, in principle, one needs to perform a full diagonalization of the hamiltonian 
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matrix in the space of (1 + 3)huj excitations. Apart from the technical difficulties, 
this solution has the disadvantage that most of the eigenstates \n) produced will 
couple only weakly to the initial and final states. 

To face the above problem, we have adopted the BAGEL approach of Skouras and 
Miither^*^) which is suitable for selecting specific eigenstates of a hamiltonian matrix 
of large dimension. We describe below the manner in which this approach was applied 
to the current problem. 

We consider the general case where the states \i) and |/) are different. Starting 
from these we can generate four states, to be denoted by {(pf^) {i = 1,..,4), by 
requiring them to satisfy the following equations: 

(0fVf)=5.., ^,J = h...A , (18) 





= ctii < 










= a2i|< 










= ctsil' 






#) 


E''\f) 


= (341 < 






i>f) + a,MT) 



From the above discussion regarding the operator E^^ it follows that the first two 
states 101°^) and 102°'') spurious. On the other hand, the other two states 103°'') 
and 104°^) have the following properties: i) they are orthogonal to the spurious states 
generated by the isoscalar El operator and ii) they are the only physical states in 
the (1 + ?))huj space which are connected by the isovector El operator to both the 
initial and final states. In the lowest order approximation of the BAGEL scheme 
the four vectors \(t)f^) {i = l,--,4) define a suitable basis for the diagonalization 
of the hamiltonian. Because of the construction of this operator, the results of the 
diagonalization retain the separation of spurious and non-spurious states. Thus the 
diagonalization provides the two physical eigenstates \n) , those with zero isoscalar 
El strength, which are to be used in (||) and (|^). 

A better approximation to the above procedure can be obtained^^) by treating 
l^l^"*) — ■■■1^) initial vectors in the Lanczos diagonalization method^^). Thus 
four more vectors, to be denoted by \(t)'P) {i = 1, ••, 4) can be constructed by requir- 
ing them to satisfy the equations: 

(0«|0f) = ^., , (0rVf) = O ^,J = 1A , (20) 
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H\<^)=T.bfl'\<^) + Y.4?\<Pr), ^ = WA . (21) 



As is evident from (pO|) the four new vectors \(t)i^) (z = 1,4) carry no El strength. 
However, as (|2T|) imphes, these are the only four vectors that are connected with \(t)f^) 
through the hamiltonian. Thus diagonahzation of H in the eight- dimensional space 
of \(t)f^) and |0f ^) (z = 1, 4) will result in a better determination of the eigenvalues 
En in (^ and (|^ while the El strength will be shared among the eigenvectors. 
The new diagonahzation again preserves the separation of spurious and non-spurious 
states while the eigenvectors corresponding to physical states are again characterized 
by having no isoscalar El strength. 

The above procedure can be continued in the manner indicated by ( pO|) and (|2T1) by 
considering additional iterations where at each iteration four additional basis vectors 
are introduced. However these new vectors interact directly only with the vectors of 
the previous iteration -in the basis of |</'j-'^'') the hamiltonian matrix is tridiagonal^^)- 
and thus only a small number of iterations is required to generate the more important 
states \n) in (^ and (^. Some examples of the convergence of expression (|^ with 
respect to the number of iterations are discussed in section 3. 

2.3 Two-body interaction and oscillator parameter 

In this section we discuss the determination of the two-body interaction G that ap- 
pears in eq. (p!oD and also the manner in which we selected the value of the oscillator 
parameter h = (h/mcuY^'^ which is used to define the basis of single-particle states as 
presented in ([T^). 

The matrix elements of G have been determined by solving the Bethe-Goldstone 
equation 

directly in the basis of harmonic oscillator states^°). For the starting energy Eg a 
constant value of -30 MeV has been adopted while the Pauli operator Q was defined 
to exclude any intermediate two-particle configuration with one nucleon in a Os or 
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both nucleons in a valence configuration, which is taken into account in our shell- 
model calculation. In eq. (^) the energy of the intermediate states is determined 
considering only the effects of the kinetic energy operator t. 

For the bare NN interaction V in (|2^ ) we have considered two versions of the Bonn 
OBE potentiaF^). The parameters of these potentials have been adjusted to fit the 
NN scattering phase shifts by solving the Thompson equation. The two potentials 
are denoted as A and C in table A. 2 of and they mainly differ in the strength 
of the tensor component. The potential C contains a moderate tensor component, 
which yields a D-state probability for the deuteron of Pd = 5.53%, while potential 
A contains an even weaker tensor component which results in Pd = 4.47%. In the 
following we shall denote the G-matrices corresponding to potentials A and C by 
and G*", respectively. 

The above potentials A and C have been used in Dirac-Brueckner-Hartree-Fock 
calculations for light nuclei^^). Such calculations have shown that the Dirac spinors 
for nucleons in nuclear medium are substantially different from those of free nucleons. 
The ratio of large to small components for the Dirac spinors in the medium may be 
described in terms of an effective mass m*. It has been shown^^'^^) that the value of 
m* = 630MeV is a reasonable choice for light nuclei. In our calculation we use this 
value of m* as well as the value m* = 938MeV of the free nucleon (which means that 
a change of the Dirac spinors in the medium is ignored) and we shall distinguish the 
corresponding G-matrices by Gm and G, respectively. 

The other parameter that enters our calculation is the oscillator parameter b. 
This is treated as a variational parameter and we adopt the value of b for which the 
binding energy of each nucleus is a minimum. This procedure has been repeated for 
each G-matrix considered in the calculation. 

3 Results of the calculation 

As outlined in sect. 2.3, in our calculation we have considered four types of two-body 
matrix elements. These correspond to using as V in eq. (|2^) versions A and C of 
the Bonn potentiaF^) and values of 938 MeV and 630 MeV for the parameter m*. 
In the following we shall distinguish the results corresponding to these four sets by 
G^, G"", G^ and G^. It should be noticed that the use of an effective mass implies 
that the structure of Dirac spinors, which are used to evaluate the matrix elements of 
the OBE potential, are modified. However, once the determination of the two-body 
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matrix elements is accomplished, the remaining part of the structure calculation is 
performed in a non-relativistic manner. This implies that we do not consider the 
effects of m* on the kinetic energy operator in (p!0|). 

All the G-matrices described above were determined for the values of 1.6, 1.8, 
2.0 and 2.1 fm for the oscillator parameter b. Thus, altogether, 16 sets of two-body 
matrix elements were determined and the properties of the low- lying states of ^Li 
and ^Li were calculated for all these interactions. Since a calculation of ^^Li in the 
large model space under consideration requires a large amount of computer time a 
more restrictive selection has been made for this isotope. 

Table 1 shows the dependence of the binding energy on the interaction and the 
oscillator parameter. The experimental values for this quantity are^^) -31.99, -39.25 
and -45.54 MeV for ^Li, ''Li and ^^Li, respectively. As seen in this table, all interactions 
considered in the calculation greatly underestimate the binding energies of the Li 
isotopes. In general more binding energy could be obtained by increasing the value 
for the starting energy Eg in the Bethe-Goldstone equation (|2^). This would imply, 
that one would treat the starting energy as a free parameter. It has been our aim, 
however, to show how far one can reproduce the empirical data without adjustment of 
parameter. Therefore we have kept the value of Eg at -30 MeV, which is a reasonable 
average of values determined in a self-consistent scheme. 

At first sight there seems to be a discrepancy between the severe underbinding 
obtained in the present investigation of Li isotopes and the very encouraging results 
on the binding energy of ^^O obtained by Skouras et al also employing OBE 
potentials and accounting for 2 hu configurations. One has to keep in mind, however, 
that the OBE potentials employed in the present investigation are different from 
those used in reference 38 although both are defined by Machleidt^^). As we have 
already mentioned above the OBE potentials used here have been obtained by solving 
the Thomson equation, while those used before^^) employ the Blankenbecler Sugar 
equation for the NN scattering equation. Furthermore the results presented here are 
obtained from a diagonalisation in a limited space of configurations. This emplies that 
the calculated energies contain effects of unlinked diagrams, which yield repulsion, as 
compared to the effective operator approach discussed in ref. 38. 

The differences obtained for the various interactions can be understood as fol- 
lows: (i) The OBE potential with a weaker tensor component (A) yields more bind- 
ing energy than the one with a stronger tensor component (C). This is a general 
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feature of BHF calculations for closed shell nuclei, employing phase-shift equivalent 
potentials'^), which obviously is valid also for binding energies obtained in shell-model 
calculations for open shell nuclei, (ii) The modification of the Dirac spinors in the 
medium contained in Gm reduces the calculated binding energy. This is also a feature 
already present in DBHF calculations and which can be understood as a reduction 
of the attractive components due the exchange of a scalar meson (a), (iii) A smaller 
calculated binding energy is accompanied by a larger value of the oscillator parameter 
b for which the minimum is obtained. In the case of ^Li and ''Li one observes that 
the minimum value for occurs for fe = 1.8 fm although an almost equal value is 
obtained for 6 = 1.6 fm. Thus a proper variational calculation using interaction 
most certainly will find the minimum between the above two values of the oscillator 
parameter. A different behavior is observed with the other three interactions where, 
as the results of table 1 indicate, the minima are shifted to larger oscillator values. 
This is particularly evident in the case of the G^ results. This behavior is similar to 
features of DBHF calculations for ^^O, in which one finds that a larger tensor force 
and the inclusion of Dirac effects increase the calculated radii. 

Considering the results of table 1 as those of a restricted variational calculation 
we obtain a natural selection for the oscillator parameter to be considered in the 
rest of the calculation. Thus for investigations on ^Li and ''Li employing G^, G'~^ 
or G^ interactions we adopt 6 = 1.8 fm, while for the interaction we consider 
h = 2.0 fm. The corresponding values for ^^Li are b = 2.0 fm and b = 2.1 fm for 
G^ and interactions, respectively. 

Fig. 1 shows the experimental and theoretical spectra of the low-lying positive 
parity states of ^Li. The theoretical predictions on the electromagnetic properties 
of this nucleus, determined for all interactions discussed above, are compared to the 
experimental data'^^) in table 2. 

As may be seen in fig. 1, the excitation energies of the ^Li states with isospin 
T = are rather insensitive on the interaction used and show a good agreement with 
the experimental data, whereas the position of the states with isospin T = 1 relative 
to the T = states changes quite drastically with the choice of the interaction. This 
behavior is particularly evident in the case of the J = 0, T = 1 state where one 
obtains excitation energies which differ by almost 2 MeV. The excitation energy of 
this state is underestimated using and the agreement with the experiment becomes 
worse using G"", G^ and G^. 
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The four interactions produce quite similar results for most of the electromagnetic 
observables displayed in table 2. These results are also in reasonable agreement with 
the experimental data, bearing in mind that in our calculation we consider bare 
electromagnetic operators. One, however, observes that the calculation predicts 
Tms and B{E2) values which are considerably larger than those obtained with the 
other interactions. This behavior can only partly be attributed to the larger value of 
the oscillator parameters b used in the calculation. Therefore we conclude that 
the increase in the calculated radii as we go from interaction G"^ to describes the 
effect already observed in DBHF calculations for closed shell nuclei (see discussion 
above). 

Fig. 2 shows the experimental and theoretical spectra of the low-lying negative 
parity states of '^Li. Unlike the case of ^Li, discussed above, the four interactions 
considered in the calculation produce spectra which are quite similar to each other 
and also with the experimental one. It seems that the differences observed in the 
T = 1 spectrum as compared to the T = states of ^Li arc getting less important 
for systems with more valence nucleons. In particular, the spectra obtained with the 
G-matrices derived from potential A , independent of the choice of the parameter m*, 
are in close agreement with the data. 

To calculate the dipole polarizability for ^Li one needs to calculate the unnatural 
parity states which are connected by the El operator to both the ground and first 
excited states of this nucleus. These states have spins 1/2+, 3/2+ and 5/2+ while 
their isospin can be either 1/2 or 3/2. As discussed in sect. 2.2, these states are 
determined in the space of (1 + 3)huj excitations using the BAGEL approach. This 
approach selectively determines the states which have strong El coupling to the 
ground and first excited states. The results of this calculation are shown in fig. 3 for 
the T = 1/2 states and fig. 4 for the T = 3/2. The results shown in figs. 3,4 show 
both the energy position as well as the B[E1, n ^ i) strength from various unnatural 
parity states n to the ground state {i = 3/2~) or first excited state (i = 1/2^) and 
correspond to the interaction. As discussed below, for this interaction one obtains 
the largest, in magnitude, values of the Tjj tensors. 

As figs. 3 and 4 show, the El strength is distributed over a large number of states 
in the energy range of 10 to about 70 MeV. One may also observe a spin dependence 
in the distribution. For example, in the case of the 1/2+ states about 40% of the 
strength is concentrated in the lowest state at about 9.5 MeV, while in the case of 
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the 3/2+ states the strength is distributed over many states and one has to extend 
to about 25 MeV to exhaust 50% of the total strength. Moreover, one observes that 
there is similarity between T = 1/2 and T = 3/2 distributions although the latter is 
shifted by about 3 MeV. 

In table 3 we list the theoretical predictions on the electromagnetic properties 
of ^Li for all interactions considered in the calculation. Table 3 also includes the 
predictions of the calculation on the polarizability terms P, th and ri2. One should 
remark at this point that with the existing data on ^Li it is not possible to obtain 
an estimate for P ^) and thus the theoretical predictions for this quantity cannot be 
yet be compared with experiment. One should also observe that the sign of tu is 
arbitrary since, as is evident from eq. (^, it depends on the relative sign of the initial 
and final states which cannot be determined uniquely. The values of T12 listed in table 
3 have been determined according to the convention of Barker and Woods^^), i.e by 
assuming that (1||£'2||2) is positive. 

As may be seen in table 3, the four calculations produce very similar values for 
the magnetic moment of the ground state as well as for the B{M1) corresponding 
to the deexcitation of the first excited state. On the other hand, one observes a 
dependence of the r^s, the E2 matrix elements and the dipole polarizability terms 
on the interaction employed. This behavior is more pronounced in the results, 
but one should remember that these were obtained with a larger b value. Generally, 
the results are in closer agreement with experiment than those obtained with the 
other three interactions. 

The most sensitive quantities, of those listed in table 3, are the polarizability 
tensors rn, Tu and the monopole term P since their determination involves both 
natural and unnatural parity states. The energy spectra of the latter states have 
been determined in the BAGEL approach. Table 4 presents the dependence of the 
polarizability quantities on the number of BAGEL iterations. The results shown in 
table 4 correspond to the interaction but a similar behavior was observed for the 
other interactions. As this table shows, there is a very fast convergence of Th, T12 
and P with respect to the number of iterations. In the actual calculation of these 
quantities we considered 25 iterations but one observes from table 4 that after the 
fifth iteration the change of results is insignificant. 

As eq. indicates, the contribution of the various J, T unnatural parity states 
adds up coherently in the case of the monopole polarizability P. On the other hand. 
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as is evident from (|^, there is a dependence of the Xj/ tensors both on the spin as well 
as on the El strength of the contributing unnatural parity states. This dependence 
is examined in table 5 for both and interactions. As table 5 shows there is a 
coherent contribution from all J, T states to T12 with the J = 1/2+ states contributing 
more. In contrast there are cancellations in the various contributions to tu. It is 
interesting to observe that for a Jj = 3/2 ground state, as is the case of ''Li and ^^Li, 
Til is given by 

^(_l).„+3/2^(n^| 2^) ^(^''"^/^ (23) 

with the Racah coefficients being positive for all possible J„ values. Thus the cancel- 
lations in Til are entirely due to the phase factor in (^). 

We conclude our study of ''Li by examining the effects of configuration space on 
the properties of this nucleus. To study these effects we made additional calculations 
in the space of Ohu configurations for the natural parity states and Ihuj for the 
unnatural parity ones. In fig. 5 we compare the low-lying energy spectrum obtained 
in the restricted space with that of the extended space while table 6 lists the remaining 
properties of the natural parity states. As the results shown in both fig. 5 and table 
6 indicate, there is a drastic improvement in the theoretical predictions when one 
extends the model space. One, of course, should note that the restricted space results 
could also be improved by a renormalization of the operators^^). 

The polarizability quantities T12, tu and P have also been calculated for different 
choices of model spaces and the results are shown in table 7. In one calculation, termed 
"small" in the table, the natural parity states have been calculated in the Ohu space 
and the unnatural in the Ihu space. In another calculation, termed "medium" in 
table 7, the space of the natural space is extended to (0 + 2)huj excitations, while the 
unnatural parity states are again determined in the Ihuj space. Finally, "large" in 
table 7 denotes the results obtained in the complete space used in this calculation. 

As may be seen in table 7, the monopole polarizability P obtains its maximum 
value in the small space. On the other hand, the values obtained in this space for the 
two quadrupole tensors are about 50% in magnitude of those obtained in the complete 
space. The smallest values in magnitude are obtained for all three quantities in the 
medium space calculation. The reason for this behavior is twofold: i) the coupling 
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between 2huj and Ihuo configurations is weak and ii) there is a considerable increase 
in the energy denominators in eqs @ and (|^ as can be deduced from the increase of 
binding energies hsted in table 6. 

The results in both tables 6 and 7 clearly suggest the importance of high-lying 
configurations which affect both the properties of the low-lying states, as well as the 
polarizability effects. Extrapolating this behavior one expects a further improvement 
in the shell-model results if even higher configurations, like for the natural parity 
states and bTioj for the unnatural parity ones could be included. Such an enlargement 
of the model space could most probably make the predictions of the shell-model 
similar to those of the cluster model. 

The nucleus ^^Li is known^^) to be a loosely bound system with a very small 
two-neutron separation energy. Therefore, the need for considering a very large shell- 
model space would be more pronounced for this nucleus than for ''Li. Hence, we do not 
expect our parameter free calculation to account for all the properties of ^^Li despite 
the fact that we employ a space in which all configurations up to "ihuo excitation are 
included. 

From the available experimental data on ^^Li, one knows that this nucleus has 
a very large r^s value of 3.16 ± 0.11 fm ^'^) and in addition there is evidence that 
the first excited state at 1.2 MeV has positive parity^^). From the ^Li investigation, 
described above, we know that such effects are best described in our model if one 
uses the and interactions. Therefore, the results to be discussed below were 
obtained with the use of these two interactions. 

Fig. 6 shows the predictions of our calculation on the low-energy spectrum of ^^Li. 
As this figure shows the calculation predicts that the first excited state of ^^Li is a 
1/2— state followed by a series of positive parity states the lowest of which, a 3/2+, 
appears at about 4 MeV. As discussed above, one expects that the results shown in 
fig. 6 could change considerably by the inclusion of higher configurations in the model 
space. 

Fig. 7 shows the energy position and the distribution of El transition strengths 
to the ground state for the T = 5/2 positive parity states of ^^Li. The corresponding 
information for the T = 7/2 states is shown in fig. 8. The results shown in figs. 
7 and 8 have been obtained using the interaction which, as evidenced from fig. 
6, produces lower excitation energies for the unnatural parity states. A comparison 
of the distributions shown in figs. 3 and 7 shows a significant difference between 
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^Li and ^^Li. In the latter, particularly in the 3/2"'" case, one observes a low-energy 
component in the distribution. Specifically, the lowest 3/2"'" state, predicted to be at 
3.86 MeV, carries about 4% of the total El strength. This feature could be interpreted 
to correspond to the soft dipole mode speculated for this nucleus^^"^^). 



Finally, in table 8 we summarize the predictions of our calculation regarding the 
ground state properties of ^^Li. As this table shows the calculation accounts satisfac- 
torily for the Tms and the magnetic and quadrupole moments of this nucleus. This is 
particularly true for the Vms value obtained with the interaction. It is interesting 
also to note in table 8 that the calculation predicts quite larger P values than those 
obtained for ^Li with the same interactions. This should be attributed to a) that 
a larger b value is used in the ^^Li calculation and b) that the excitation energy of 
the positive parity states is greatly reduced. On the other hand the ^^Li calcula- 
tions predict smaller rn values than those obtained for ^Li. This behavior should 
be attributed to the strong cancellation among the contributions of the various J, T 
positive parity states, an effect discussed above in connection with ''Li. It would be 
of considerable interest if values of P and th were obtained in future experiments on 
^^Li, since these will provide a useful test of our calculation. 



Results of shell- model calculations for the isotopes *^Li. 'Li and ^ILi are presented, 
which consider configurations within various major shells. Realistic hamiltonians 
are considered, which contain the kinetic energy and a NN interaction derived from 
modern OBE potentials'^). The effects of NN short-range correlations are taken into 
account by solving the Bethe-Goldstone equation for these potentials, considering 
a Pauli operator which is consistent with the shell-model configurations taken into 
account. No further renormalization of the hamiltonian and the operators for the 
electromagnetic transitions has been made since it is the aim to account for those 
long-range correlations by a sufficiently large shell-model space. The effects due to 
spurious center of mass motion are considered by adding an appropriate harmonic 
oscillator potential for the center of mass coordinate. The main conclusions can be 
summarized as follows: 

• The bulk properties (binding energies, radii) calculated for the open shell nuclei 
show a similar dependence on the OBE interaction used as it has been observed 
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in DBHF calculations for closed shell nuclei: NN interactions with a stronger 
tensor component yield less binding energy as a phase-shift equivalent potential 
with a weaker tensor force; the modification of the Dirac spinors for the nucleons 
in the medium reduces the calculated binding energy; a smaller binding energy 
is correlated with a larger value for the radius. 

The calculated excitation spectra are weakly depending on the NN interaction. 
Only in the case of ^Li a strong dependence of the energies for the states with 
isospin T = 1 relative to those with T = is observed. A good agreement with 
the experimental data is achieved if a large model space is considered. 

Also the calculated electromagnetic properties show a good agreement with the 
empirical data, keeping in mind that the present calculation does not contain 
any adjustable parameter. Of course it is evident that our predicitons are not 
as good as those obtained in phenomenological studies employing parametrized 
effective hamiltonians and electromagnetic operators. In particular our results 
for the radii show poor agreement with the empirical values. This is a well 
known problem of microscopic studies of light nuclei, employing modern OBE 
potentials 1^ 

The results of the present investigation clearly indicate (see table 7) that the 
polarizability tensors tu and T12 for ^Li depend strongly on the model space. 
The fact that the present calculation produced values for these quantities which 
are much closer to the experimental estimates than those obtained in previous 
shell-model studies^^"^^) should be, therefore, mainly attributed to the use of 
a larger space. Thus it appears that to improve further the agreement with 
the experimental data one needs to go beyond those model spaces considered 
in the present approach. Such an expansion of the model space is currently 
very difficult to attempt due to the exceedingly large number of shell-model 
configurations involved. 

The evaluation of the dipole polarizability requires a detailed information on 
the El excitation of states which appear in a wide energy interval (see fig- 
ures 3,4,7 and 8). It was demonstrated that the BAGEL approach provides a 
powerful tool to account for such contributions in a very efficient and reliable 
way. This indicates that the same approximation could be very useful for mi- 
croscopic studies of double (3 decay or double charge exchange processes, where 
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similar summations over intermediate states occur. It is important to repeat at 
this point that the BAGEL method corresponds to exact diagonahzation pro- 
vided the number of iterations is chosen to match the dimension of the energy 
matrix^^). However, the advantage of the method is that one can terminate the 
calculation once convergence of the quantity under consideration has been ob- 
tained. It was found in the present calculation (see table 4) that a small number 
of iteration is required to obtain convergence for the polar izability tensors. On 
the other hand, a larger number of iterations will most probably be required in 
cases of non-collective excitations in which many states need to be considered 
each participating with a small transition strength. 

These investigations have been supported by the Commission of the European 
Communities, the British Council in Athens and the Bundesministerium fiir Forschung 
und Technologie (Germany, Project 06 Tii 736). This support is gratefully acknowl- 
edged. 
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Table 1 

Dependence of the binding energies of the Li isotopes on the oscillator parameter b 



b{fni) 


1.(3 


1.8 


2.0 


2.1 






^Li 










-ZU.OO 


1 Q ZLQ 






-14.46 


-15.85 


-15.43 






-15.04 


-16.64 


-16.49 


-15.93 




-8.79 


-11.48 


-12.17 


-12.00 






^Li 






qA 


-26.47 


-26.59 


-24.63 




qC 


-18.89 


-20.51 


-19.75 




r<A 


-19.54 


-21.56 


-21.24 


-20.45 




-11.77 


-15.26 


-16.03 


-15.74 






iiLi 






r^A 






-22.03 


-21.68 








-13.87 


-14.32 
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Table 2 

Electromagnetic properties of ^Li determined with various G-matrices 



Quaiitit}' 




Units 


Expeiiniont") 




G'' 


r<A 


r<C 






fm 


2.09 ±0.02 


2.34 


2.39 


2.40 


2.64 








0.8220 


0.8530 


0.8453 


0.8540 


0.8524 


Qs 




efm^ 


-0.083 


-0.302 


-0.407 


-0.099 


-0.365 


B{E2;3,0^ 


1,0) 


WU 


16.5 ± 1.3 


4.96 


5.49 


5.57 


8.35 


S(M1;0,1 ^ 


1,0) 


wu 


8.62 ±0.18 


8.64 


8.52 


8.84 


8.72 


B{E2;2,0^ 


1,0) 


WU 


6.8 ±3.5 


4.06 


4.24 


5.26 


7.04 


B{M1;2, 1 ^ 


1,0) 


WUxlO-2 


8.35 ± 1.5 


0.04 


0.09 


2.2 


0.06 



) Ref. 43). 



Table 3 

Electromagnetic properties of ^Li determined with various G-matrices 



Quantity 


Units 


Experiment") 


QA 












fm 


2.23 ±0.02 


2.41 


2.45 


2.46 


2.71 


n 






3.2564 


3.0789 


3.0918 


3.1270 


3.1339 


B{Ml-l/2 - 


-3/2) 


wu 


2.75 ±0.14 


2.39 


2.36 


2.39 


2.38 


B{E2; 1/2 - 


-3/2) 


WU 


19.7 ± 1.2 


7.72 


8.74 


9.25 


13.7 


B{E2;7/2- 


-3/2) 


wu 


4.3 


3.69 


4.01 


5.65 


6.15 


Qs 




efm^ 


-4.00 ±0.06 


-2.54 


-2.70 


-2.81 


-3.39 


B{E2;3/2 ^ 


► 1/2) t 




7.27 ±0.12 


3.14 


3.55 


3.76 


5.56 


P 




fw? 




0.095 


0.109 


0.110 


0.164 


Til 




fm^ 


-0.12 ±0.07 


-0.050 


-0.056 


-0.057 


-0.084 


~12 






-0.118 ±0.012 


-0.019 


-0.055 


-0.053 


-0.080 



") The data for the last five quantities are from ^). The rest are taken from ^3). 
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Table 4 

Dependence of the dipole polarisabitity on the number of BAGEL iterations 



Iteration 


' 12 


' 11 


P 


1 


-0.06377 


-0.06662 


0.13664 


2 


-0.07759 


-0.08138 


0.16131 


3 


-0.07890 


-0.08368 


0.16386 


4 


-0.08067 


-0.08490 


0.16459 


5 


-0.07968 


-0.08325 


0.16406 


6 


-0.07997 


-0.08367 


0.16417 


7 


-0.08004 


-0.08376 


0.16416 


8 


-0.08005 


-0.08377 


0.16415 


9 


-0.08005 


-0.08376 


0.16414 


Final 


-0.08005 


-0.08376 


0.16413 



Table 5 

Contribution of different J, T states to Tj/ 



T 


J'^ 


Tl2 


Til 


G 

Tl2 


C 

m 

Til 


1/2 


1/2+ 


-0.03465 - 


■0.05039 


-0.05220 


-0.07311 




3/2+ 


-0.00755 


0.02759 


-0.01113 


0.04062 




5/2+ 




■0.02091 




-0.03072 


3/2 


1/2+ 


-0.01038 - 


■0.02802 


-0.01630 


-0.04351 




3/2+ 


-0.00001 


0.02958 


-0.00042 


0.04531 




5/2+ 




-0.01189 




-0.02235 


Total 


-0.05259 - 


■0.05702 


-0.08005 


-0.08376 
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Table 6 

Effects of configuration space on tlie electromagnetic properties of the low-lying states of ^Li 



Quantity 


Units 


Experiment") 






















(0 + 2)hu 




(0 


+ 2)hu; 


B.E 




MeV 


39.25 


14.11 


21.24 


9.59 




16.03 






fm 


2.23 ±0.02 


2.50 


2.46 


2.78 




2.71 








3.2564 


3.1431 


3.1270 


3.1339 




3.1375 


5(M1; 1/2 - 


-3/2) 


WU 


2.75 ±0.14 


2.45 


2.39 


2.45 




2.38 


B{E2;l/2 - 


-3/2) 


wu 


19.7 ± 1.2 


4.17 


9.25 


6.67 




13.7 


B{E2;7/2 - 


-3/2) 


WU 


4.3 


1.64 


4.06 


2.61 




6.15 


Qs 




efrn? 


-4.00 ±0.06 


-2.00 


-2.81 


-2.46 




-3.39 




■ 1/2) T 




7.27 ±0.12 


1.70 


3.76 


2.71 




5.56 



") The data for the last two quantities are from ^). The rest are taken from ^^). 



Table 7 

Effects of configuration space on the dipole polarizability of ^Li 



Space 


Ti2 ifm^) 


Til ifm^) 


P ifm^) 


Small") 


-0.0431 


-0.0463 


0.197 


Medium*") 


-0.0409 


-0.0432 


0.095 


Large*^) 


-0.0801 


-0.0838 


0.164 


Experiment*^) 


-0.148 


-0.12 





") Ohcu space for the negative parity states and Ihuj for the positive parity states. 
'') (0 + 2)^a; space for the negative parity states and Ihcu for the positive parity states. 

(0 + 2)hu! space for the negative parity states and (1 ± 3)hu! for the positive parity 
states. 
'^) Ref.^). 
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Table 8 

Properties of ^^Li 



Quantity Expor.") G;^ G^,; 

rmsifm) 3.16 ±0.11 2.81 2.97 

fi (fi^) 3.6673 ± 0.0025 3.6908 3.6890 
Qs (efm^) -3.12 ± 0.45 -3.48 -3.92 

Til (/m^) -0.0149 -0.0624 

P (fm^) 0.238 0.375 



") The value for is from . For the other quantities from '^^) . 
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Figure Captions 

Fig. 1: Low-lying spectra of ^Li calculated with various G-matrices (see text) in the 
space of (0 + 2)lyjj excitations. All states have positive parity. Each level is 
labelled by J, T. The experimental information is from Ref. ^^). 

Fig. 2: Low-lying spectra of ^Li calculated with various G-matrices (see text) in the 
space of (0 + 2)/ia; excitations. All states have negative parity and T = 1/2. 
Each level is labelled by 2 J. The experimental information is from Ref. '^'^). 

Fig. 3: Distribution of the El strength among the T = 1/2 positive parity states 
of "^Li. Only transitions leading to the ground state 3/2^ and to the first 
excited 1/2^ state are shown. The accumulated strength is normalized to 
1 for the complete strength in each channel. 

Fig. 4: Distribution of the E\ strength among the T = 3/2 positive parity states 
of ''Li. Only transitions leading to the ground state 3/2~ and to the first 
excited l/2~ state are shown. 

Fig. 5: Low-lying spectra of ^Li calculated with interactions Gf^ and (see text). 

The spectra labelled (a) have been determined in the QTiuj space, while 
those labelled (b) in the (0 -|- 2)hu space. All states have negative parity 
and T — 1/2. Each level is labelled by 2 J. The experimental information 
is from Ref. ^^). 

Fig. 6: Low-lying spectra of ^^Li calculated with interactions Gj^ and (see text). 
All states have T = 5/2. Each level is labelled by 27^^. 

Fig. 7: Distribution of the El strength among the T — b/2 positive parity states 
of ^^Li. Only transitions leading to the ground state are shown. 

Fig. 8: Distribution of the El strength among the 7=7/2 positive parity states 
of ^^Li. Only transitions leading to the ground state are shown. 
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